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ABSTRACT 

We develop a numerical scheme for obtaining the r-mode oscillations of rapidly 
rotating relativistic stars. In the present scheme, we neglect all metric perturbations 
and only take account of the dynamics of the fluid in the background spacetime of the 
unperturbed star (the relativistic Cowling approximation). We also assume the star is 
barotropic, i.e. neutrally stable against convection under the assumption of adiabatic 
oscillations. Our numerical scheme is based on the Yoshida-Eriguchi formulation for 
the analysis of the general relativistic f-mode oscillations in the Cowling approximation 
and a general relativistic generalization of the Karino-Yoshida-Yoshida-Eriguchi's nu- 
merical scheme for obtaining oscillations of rapidly rotating Newtonian stars. By this 
new numerical scheme, frequencies of the r-mode oscillations are obtained as functions 
of the ratio of the rotational energy to the absolute value of the gravitational energy 
T/|W] along sequences of poly tropic equilibrium stars whose ratio of the pressure to 
the total energy density at the center of the star and the polytropic index are kept 
constant. It is found that the dimensionlcss oscillation frequency a/Cl is a linearly 
decreasing function of T/|W|, where a and tt are the oscillation frequency and the 
angular velocity of the star measured in an inertial frame at spatial infinity. We also 
find that oscillation frequencies of the r-modes are highly dependent on the relativistic 
factor M/R of the star as already found in previous studies in which the slow rotation 
approximation has been used. Here M and R denote the mass and radius of the star, 
respectively. 

Key words: relativity - stars: neutron - stars: oscillations - stars: rotation 



1 INTRODUCTION 

The inertial modes as well as the /-modes are basic oscil- 
lation modes of rotating stars in the sense that they exist 
even in incompressible stars (Bryan 1889; Lindblom & Ipser 
1999). The inertial mode oscillations are mainly restored 
by the Coriolis force due to stellar rotation and their os- 
cillation frequencies are therefore comparable to the stellar 
rotation frequencies (Papaloizou & Pringle 1978; Provost, 
Berthomieu & Rocca 1981; Saio 1982; Unno et al. 1989). 
The r-mode oscillations belong to a sub-class of the iner- 
tial mode oscillations and are characterized by small radial 
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displacement vectors of the eigenfunctions. They are unsta- 
ble against gravitational radiation reactions for all rotating 
stars if the stars are inviscid (r-mode instability: Andersson 
1998; Friedman & Morsink 1998). The r-mode instability is 
attributed to the so-called Chandrasekhar-Friedman-Schutz 
(CFS) mechanism, which makes various modes of oscilla- 
tions of rotating stars unstable (Chandrasekhar 1970; Fried- 
man & Schutz 1978). The r-mode oscillations excited in neu- 
tron stars due to the CFS mechanism may slow the rotation 
of newly born neutron stars as well as cold old neutron stars 
in low mass x-ray binary (LMXB) systems. Furthermore, 
there have been suggestions that the gravitational radiation 
emitted due to the excitation of the r-mode oscillations will 
be observed by the Laser Interferometric Gravitational Wave 
Observatory (LIGO) and other detectors (Lindblom, Owen 
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& Morsink 1998; Andersson, Kokkotas & Schutz 1999; Owen 
et al. 1998; Bildsten 1998; Andersson, Kokkotas & Ster- 
gioulas 1999; See also the review by Andersson and Kokkotas 
2001), though the non-linear saturation of the mode may in- 
hibit the mechanism to work in newly born stars (Arras et 
al. 2003). 

So far, most studies of the r-mode oscillations of neu- 
tron stars have been carried out within the slow rotation ap- 
proximation and/or the weak gravitational field approxima- 
tion (Newtonian dynamics). However, analyses under these 
approximations are not sufficient to investigate the r-mode 
oscillations of neutron stars because of the following reasons. 
First, the effect of general relativity on the r-mode oscilla- 
tions of neutron stars would be significant because the rel- 
ativistic factor of neutron stars, M/R, could be as large as 
M/R ~ 0.2, where M and R denote the mass and radius 
of the star, respectively. Since the effect of dragging of iner- 
tial frame on the angular velocity of the fluid measured by 
the locally non-rotating observer is in this order, the iner- 
tial mode should be also modified to this order (see Ea. Qlfil 
and the discussion there). Second, from time to time, the 
r-mode instability of nascent neutron stars has been investi- 
gated because this is one of the most important issues of the 
r-mode instability related to neutron stars. The spin rates 
of neutron stars at their birth are believed to be near the 
mass-shedding limit where the mass of the star begins to 
shed from the stellar surface of the equator. If nascent neu- 
tron stars are considered, therefore, the effect of deformation 
of the star due to rapid stellar rotation or of the centrifu- 
gal force on the r-mode oscillations should not be neglected. 
In short, both effects of general relativity and rapid rota- 
tion have to be taken into account in studies of the r-mode 
oscillations of neutron stars. 

The r-mode oscillations of rapidly rotating compress- 
ible stars have been obtained for the first time by Yoshida 
et al. (2000) for uniformly rotating stars within the frame- 
work of Newtonian gravity. They have shown that the basic 
properties of the r-mode oscillations of slowly rotating stars 
do not change greatly and that extrapolating the results of 
the slow rotation approximation gives us good estimates for 
properties of the r-mode oscillations and the r-mode insta- 
bility even for the rapidly rotating models (see, also, Karino 
et al. 2000). The Newtonian scheme has been applied to dif- 
ferentially and rapidly rotating stars by Karino, Yoshida & 
Eriguchi (2001). 

As for general relativistic investigations of the r-mode 
oscillations, Kojima (1998) considered for the first time the 
axial parity perturbations of slowly rotating stars by as- 
suming the oscillation frequencies to be proportional to the 
stellar angular velocities (low-frequency approximation) and 
derived a single second order ordinary differential equa- 
tion (Kojima's equation) for the r-mode-like oscillations. 
He found that the relativistic r-mode oscillations are ab- 
solutely different from Newtonian ones in the sense that 
Kojima's equation allows a continuous spectrum for the 
eigenfrequency (see, also, Beyer & Kokkotas 1999). Lock- 
itch, Andersson & Friedman (2001) noticed that Kojima's 
equation is valid only if stars are non-barotropic and that, 
for barotropic stars, the polar parity perturbations have to 
be taken into account as well as the axial parity pertur- 
bations in order to obtain the r-mode-like oscillations. For 
barotropic stars, they then obtained a general relativistic 



counterpart of the r-mode oscillation in Newtonian stars. 
Their eigenf unctions are composed both of the polar and ax- 
ial parity perturbations (see, also, Lockitch 1999; Lockitch, 
Friedman & Andersson 2003). Their analysis has been re- 
cently extended to superfluid stars by Yoshida & Lee (2003) . 

For the relativistic r-mode oscillations of non-barotropic 
stars, the situation is different from that of barotropic stars. 
For non-barotropic stars, as mentioned before, if we use the 
low-frequency approximation and the slow rotation approx- 
imation, basic equations becomes Kojima's equation, which 
contains a singular point inside the star for some frequency 
range and yields a continuous spectrum and singular eigen- 
functions (Kojima 1998). The main issue, which has not 
been fully solved yet, is whether there exist the r-mode- 
like solutions characterized by discrete eigenfrequencies and 
regular eigenfunctions in non-barotropic relativistic stars or 
not. In spite of various attempts to improve our understand- 
ing of the relativistic r-mode oscillations, we have not yet 
understood the r-mode oscillations of non-barotropic rel- 
ativistic stars clearly (Kojima & Hosonuma 1999; Kojima 
& Hosonuma 2000; Lockitch, Andersson & Friedman 2001; 
Yoshida 2001; Ruoff & Kokkotas 2001; Yoshida & Futamase 
2001; Ruoff & Kokkotas 2002; Yoshida & Lee 2002; Ruoff, 
Stavridis & Kokkotas 2003). The present study focuses upon 
the r-mode oscillations of barotropic stars. We will study 
non-barotropic stars in a forthcoming paper. 

The main purpose of this paper is to improve our under- 
standing of properties of the r-mode oscillations of rapidly 
rotating relativistic stars. Particularly, we want to explore 
both the effects of general relativity and rapid rotation on 
the r-mode oscillations. Up to now, there has been no pre- 
scription for managing pulsations of rapidly rotating rela- 
tivistic stars because of difficulties due to gauge conditions, 
boundary conditions at spatial infinity and other factors, 
although there is a method to obtain zero frequency /- 
modes devised by Stergioulas & Friedman (1998). However, 
this method is not available for obtaining general oscilla- 
tion modes of non-zero frequencies. In order to avoid these 
difficulties, we employ a relativistic version of the Cowling 
approximation in which all the metric perturbations are ne- 
glected (McDermott, Van Horn & Scholl 1983; Finn 1988; 
Ipser & Lindblom 1992; Yoshida & Eriguchi 1997). As dis- 
cussed by Lindblom & Splinter (1990), Yoshida & Kojima 
(1997) and Yoshida & Eriguchi (1997), the relativistic Cowl- 
ing approximation is a good approximation for determining 
frequencies of oscillation modes for which the fluid motions 
dominate vibrations of the spacetime. Therefore, the Cowl- 
ing approximation would be suitable for analyzing the r- 
mode oscillations for which the fluid motions are dominat- 
ing. 

Making use of the relativistic Cowling approximation, 
we calculate the r-mode oscillations of rapidly rotating stars. 
In this paper, we assume oscillations to be adiabatic and 
consider only barotropic stars as mentioned before. In or- 
der to obtain the r-mode oscillations of rapidly rotating 
relativistic stars, we make use of the relativistic Cowling 
formulation for the rapidly rotating polytropes by Yoshida 
& Eriguchi (1997) and generalize to rapidly rotating gen- 
eral relativistic stars the Karino- Yoshida- Yoshida-Eriguchi 
numerical method (Karino et al. 2000) for obtaining the r- 
mode oscillations of rapidly rotating Newtonian stars. The 
generalization is straightforward. Since we employ the rel- 
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ativistic Cowling approximation, basic equations are those 
of fluid motions and very similar to those of the Newtonian 
stars (we do not need to solve the linearized Einstein equa- 
tion). In §2 we present the basic equations and the numer- 
ical method employed in this paper for studying pulsations 
of rapidly rotating relativistic stars. Note that notations are 
not the same as those of Yoshida and Eriguchi (1997). Nu- 
merical results are given in §3, and §4 is devoted to dis- 
cussion and summary. In this paper, we use units in which 
c = G = 1, where c and G denote the velocity of light and 
the gravitational constant, respectively. 



2 FORMULATION 

2.1 Equilibrium State 

We consider oscillations of axisymmetric equilibrium con- 
figurations of uniformly rotating polytropic stars without 
meridional circulations. The spacetime for such configura- 
tions can be described by the line element: 



ds 2 = -e 2v dt 2 + e 2a (dr 2 + r 2 d9 2 ) 

+ e r sin 9 (dip — u> dt) . 



(1) 



where u, a, j3, and u are functions of r and 9 only. Here 
spherical coordinates (r, 9, p>) are used. The four velocity of 
the fluid, u^ , is given by 



h" =7(£ M + fi^ M ). 



VT 



(2) 



(3) 



neglect all the metric perturbations and solve only the lin- 
earized continuity equation and the linearized Euler equa- 
tions, assuming the perturbations to be adiabatic. The dy- 
namics of the fluid can then be described in terms of the rest 
mass density perturbation, Sp, and the perturbations of the 
fluid velocity, <5u M , where Sq denotes the Eulerian perturba- 
tion of a physical quantity q. Since perturbations of axisym- 
metric stars in stationary states are considered, we assume 
that the perturbations have harmonic time-dependence and 
ip- dependence given by e - tc,t + lm f ^ wn ere m is an integer 
and a is an oscillation frequency measured by an inertial ob- 
server at spatial infinity. The pulsation equations are given 
as follows. 

For adiabatic perturbations, 



p V Nip 



(6) 



where Sp is the Euler perturbation of the pressure and we 
have assumed a barotropic relation with the same index for 
both equilibrium states and for pulsations, i.e. the adiabatic 
index, F, is taken to be F — 1 + N^ 1 . Note that due to 
the barotropic equation of state, the star is neutrally stable 
against convection. 

The perturbed continuity equation is expressed as 

u* , n . r dSu r d8u e 

_( ff - m fi),S„-_-_ 



2 I dp d ,„ 

- + -JT + 1T (2a + + i/ 
r p or or 

cot6 + -p% + lT9 {2a + p + ' 



+ {<jF - to} 6u v = , 



Su r 



Su 



(7) 



J-" 



(fi-w) 



(4) 



where i M and p>^ are, respectively, the time-like and rota- 
tional Killing vectors, and fl is the constant angular velocity 
of the star. The polytropic equation of state is defined by 

r^ 1+1/iV 

p — K p ' . 



e = p + Np, 



(5) 



where p, p, and e mean the pressure, the rest mass density, 
and the total energy density, respectively. Here, JV and K 
are the polytropic index and the polytropic constant, respec- 
tively. In the present study, numerical solutions of the equi- 
librium states are constructed with the numerical method 
devised by Komatsu, Eriguchi & Hachisu (1989). 



2.2 Pulsation Equations in the Cowling 
Approximation 

In general relativity, not only the fluid but also the spacetime 
is dynamical. We therefore have to solve linearized Einstein 
equation in order to determine exact oscillation frequencies 
of relativistic stars. When we consider oscillations of rapidly 
rotating stars within the framework of general relativity, it 
is too difficult to solve the linearized Einstein equation. No 
one has succeeded in obtaining a general method for treat- 
ing pulsations in rapidly rotating relativistic stars. In order 
to avoid this difficulty, we employ a relativistic version of 
the Cowling approximation for determining the r-mode os- 
cillation frequencies of rapidly rotating stars. We therefore 



SI — LU + V 2 U) 



(«) 



where it* , 8u r = i8u r , 8u e = iSu , and 8u v , are the t- 
component of the four velocity, the Euler perturbations of 
the r-component, ^-component, and (^-component of the four 
velocity, respectively. The perturbed equations for the r- 
component, ^-component and ^-component of the equations 
of motion are: 



+ 



1 e- 2a dSp 



1 



dp 



e + p vf dr (e + p) 2 u l dr 
+ (<7 — mQ)Su r 



{5e + Sp) 



e 2 "- 2 V sin 2 9(Q - u) ?^8u* = 
dr 

1 e- 2a dSp 1 e- 2a dp 

e _|_ p r 2 u t qq (e -\- p) 2 r 2 ^ d9 

,9 



(9) 



(Se + Sp) 



+ (cr — mQ)Su 

2f3-2a ■ 2 Qln . \" m -^r f n 

— e y sin 9(il — ui) — j--— Su ^=0 



(10) 
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(11) 
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H 



fi-wl- 



(12) 



where 5e is the Euler perturbation of the energy density. 

Physically acceptable solutions of pulsation equations 
have to satisfy boundary conditions on the rotation axis and 
on the stellar surface. The regularity conditions are imposed 
on the rotation axis, i.e. all the eigenfunctions have to vanish 
on the rotation axis if non-axisymmetric perturbations are 
considered. The surface boundary condition is given by 



ry(u — mfl) Ap = —17 (a — mQ)5p + Su* 1 



dp 



0, (13) 



where Aq means the Lagrangian perturbation of a physical 
quantity q, which is related to the Eulerian perturbation 5q 
through the relation Aq — Sq + C^q. Here £5 denotes the 
Lie derivative along the Lagrangian displacement vector £ M . 
As in the Newtonian models (see e.g. Karino et al. 
2000), we introduce a surface-fitted coordinate system as 
follows: 

e, = e , (14) 



r.(6) 



where r — r s (8) denotes the surface shape of the unper- 
turbed configuration. In the actual numerical computations, 
we make use of equidistant mesh points both in the ra- 
dial and angular directions (0 < r» < 1; < 0» < 7r/2). 
Note that it is sufficient to consider the angular range 
< 0, < 7r/2 to determine oscillation modes of rotating 
axisymmetric equilibrium stars, because eigenfunctions of 
the stars must have reflection symmetry or reflection anti- 
symmetry with respect to the equatorial plane. For the r- 
mode oscillations we consider, the eigenfunctions are reflec- 
tion anti-symmetric. In the present investigation, we take 
the number of mesh points to be mostly (r*x#„) = (101x70) 
and sometimes (r» x 0,) = (121 x 101). In order to ob- 
tain numerical solutions of the r-mode oscillations, pulsa- 
tion equations described before are transformed to a new 
surface-fitted coordinate system and discretized on the mesh 
points, and are cast into coupled non-linear algebraic equa- 
tions in terms of the oscillation frequency a and the dis- 
cretized eigenfunctions, which are solved with the Newton- 
Raphson iteration scheme. The numerical procedure we use 
is based on that of Yoshida and Eriguchi (1997) and is a 
generalization of that of Karino et al. (2000). 



3 NUMERICAL RESULTS 

Since the m = 2 r-mode is the most unstable against gravi- 
tational radiation reactions among the inertial modes (Lock- 
itch & Friedman 1999; Yoshida & Lee 2000a; Yoshida & Lee 
2000b), it would be the most important mode in the spin 
evolution of neutron stars. Therefore, in the present study, 
we concentrate on the r-mode oscillations associated with 
m — 2. There are neither overtone modes nor modes associ- 
ated with I 7^ \m\ in the non-rotation limit, where I is an an- 
gular eigenvalue of the spherical harmonic function, for the 
r-mode oscillations of barotropic stars (Provost, Berthomieu 
& Rocca 1981; Lockitch & Friedman 1999). Thus, only fun- 
damental r-modes with I — \m\, whose eigenfunctions have 
no node in the radial direction, are calculated. By decreasing 
the axis ratio r p /r e of the equilibrium star from r p /r e = 1, 



Table 1. Equilibrium sequences for which the r-modc oscillations 
have been analyzed. Here, the relativistic factors M/R shown in 
this table are evaluated in the non-rotation limit. Here M and 
R arc the gravitational mass and the circumferential radius of 
the equilibrium models. Therefore, the values of M/R along the 
equilibrium sequences are not constant but their variations are 
not large (typically a few percent for N = 0.5 and 10% for TV = 1). 
p c and £ c are the pressure and the energy density at the center 
of the star. 



Sequence 


N 


M/R 


Pc/£c 


a 


0.5 


0.1 


0.06180 


b 




0.2 


0.1744 


c 


1.0 


0.1 


0.06614 


d 




0.2 


0.2016 



Table 2. i?-modc frequency in the limit of no rotation. The re- 
sults of Yoshida and Lee (2002) in which the relativistic Cowl- 
ing approximation is used are also shown. The results of the full 
problem are obtained by solving the full equations for the r-modc 
oscillations in the slow rotation approximation. The quantities u> c 
and 0J e q are the values of 
the surface, respectively. 



Method 



N 



M/R Pc/sc Weq/H w c /Q a/Q. 



Present 0.5 0.10 0.06180 0.073 0.236 1.410 

Yoshida & Lee 0.5 0.10 1.410 

Full 0.5 0.10 1.379 

Present 0.5 0.20 0.1744 0.167 0.478 1.502 

Yoshida & Lee 0.5 0.20 1.503 

Full 0.5 0.20 1.443 

Present 1.0 0.10 0.06614 0.059 0.274 1.410 

Yoshida & Lcc 1.0 0.10 1.410 

Full 1.0 0.10 — 1.380 

Present 1.0 0.20 0.2016 0.139 0.556 1.505 

Yoshida & Lcc 1.0 0.20 1.506 

Full 1.0 0.20 1.453 



where r p and r e are the polar radius and the equatorial ra- 
dius of the star, respectively, we obtain equilibrium configu- 
rations from a nonrotating spherical star to a deformed star 
at the mass-shedding limit by the Komatsu, Eriguchi, and 
Hachisu numerical scheme (Komatsu, Eriguchi & Hachisu 
1989, KEH). Along an equilibrium sequence, the ratio of 
the pressure to the total energy density at the center of the 
star, p c /Ec, the polytropic index, N, are kept constant. As 
another rotation parameter, in addition to r p /r e , we use the 
ratio of the rotational energy (T) to the absolute value of 
the gravitational energy (W), T/\W\ (for the definition of 
T/|W|, see, Komatsu, Eriguchi & Hachisu 1989). Physical 
quantities characterizing the stars along an equilibrium se- 
quence are tabulated in Table 1. 

In order to check our numerical code, we have computed 
several r-mode oscillations for slowly rotating stars and ex- 
trapolated the r-mode frequencies to the limit of £1 — > from 
the obtained frequencies. The results in the non-rotating 
limit are compared with the solutions obtained by Yoshida 
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and Lee's numerical scheme, in which all oscillation modes 
of slowly rotating relativistic stars have been computed un- 
der the relativistic Cowling approximation (Yoshida & Lee 
2002). The results summarized in Table. 2 show that oscil- 
lation frequencies obtained by two different methods are in 
excellent agreement and that the numerical scheme devel- 
oped in the present study works quite well. In this table, u> c 
and w eq denote the values of the metric function u) at the 
center and the surface of the spherical star, respectively. 

In Table. 2, the results for the r-mode oscillations of the 
full equations including the metric perturbations are also 
shown. These results are obtained by using the formalism 
proposed by Lockitch et al. (2001) (see, also, Yoshida & Lee 
2003). Although it was very difficult to obtain the smooth 
eigenfunction for the model with N = 1.0 and M/R — 0.2 
in the full analysis, the obtained eigenfrequency for that 
model should not be far from the true value. As seen from 
this table, the relativistic Cowling approximation can give 
very good approximate values for the oscillation frequencies 
to within a few percent. 

Stergioulas & Font (2001) have extracted an r-mode- 
like characteristic frequency from their general relativis- 
tic simulations. They performed hydrodynamical simula- 
tions using the relativistic Cowling approximation where the 
spacetime geometry was fixed during the time evolution of 
the fluid. Therefore, the present results can be directly com- 
pared with their results. We have obtained a frequency of 
1.05 kHz for the m = 2 r-mode oscillation of an N = 1 rotat- 
ing polytrope similar to that of Stergioulas & Font, whose 
gravitational mass, equatorial circumferential radius, and 
spin period are M = 1.63M , R = 17.17 km, and P = 1.25 
ms, respectively. Our result is in good agreement with the 
oscillation frequency of the r-mode obtained by Stergioulas 
& Font, i.e. 1.03 kHz. 

In Figures 1 and 2, the ratio of the r-mode frequency 
to the angular velocity, a/Q, of rotating stars is shown as 
a function of T/\W\ for four equilibrium sequences listed in 
Table 1. Note that in the present study the frequency mea- 
sured by an inertial observer at spatial infinity is shown; the 
corotating frequency a c used in many papers is related to a 
by cr c = a — mil. Figures 1 and 2 show results for polytropes 
with N = 0.5 (sequence a and sequence b) and N = 1.0 (se- 
quence c and sequence d), respectively. From these figures, 
it is found that the scaled eigenfrequencies of the r-mode 
oscillations, er/fi, are decreasing functions of T/\ W\ and are 
given as nearly linear functions of T/\W\ except for config- 
urations in the mass-shedding states. Similar features have 
been found in the Newtonian r-mode oscillations of rapidly 
rotating stars (Karino, Yoshida & Eriguchi 2001). It is im- 
portant to note that this simple relation between a/fl and 
T/IM 7 ! is preserved even for highly relativistic stars. Since 
the oscillation frequencies a/il are nicely approximated by 
linear functions of T/\W\ except for the region near the 
mass-shedding limit, it is meaningful to derive approxima- 
tion formulas for a/Q versus T/|W|, which are given by 

1.41 — 1.85 t^tt for sequence a 

1.50 — 1.23 t^fj for sequence b 
1.41 — 1.95 t-jtm for sequence c 

1.51 — 1.36 t^tt for sequence d. 




(15) 



apply a least-square-fit to the scaled frequencies as functions 
of T J | W | . Comparing the coefficients appearing in these for- 
mula I15H for different sequences, we can see that the r-mode 
frequencies of relativistic stars do not strongly depend on 
the stiffness of the equation of state, but they are sensi- 
tive to the relativistic factor M/R and to the rotation rate. 
Similar strong dependence of the relativistic factor M/R 
on the oscillation frequencies of the r-mode oscillations has 
been already observed in the r-mode oscillations of slowly 
rotating relativistic stars (e.g., Kojima 1998; Lockitch, An- 
dersson & Friedman 2001; Yoshida 2001; Ruoff & Kokkotas 
2001; Yoshida & Lee 2002; Lockitch, Friedman & Ander- 
sson 2003). As discussed by Kojima (1998), frequencies of 
the r-mode oscillations in the limit of no rotation can be 
approximated by a certain average of the local oscillation 
frequency of the r-mode oscillation, namely 



a 



& {" 



+ 2)(M-l)+2 



-1 



(16) 



In equation 1168 . the relativistic effect appears only in the 
second term in the braces, which is always positive for non- 
pathologically rotating stars. Therefore, the frequency of the 
r-mode oscillation tends to increase as the relativistic fac- 
tor M/R is increased because the dragging of inertial frame 
becomes significant for highly relativistic stars as seen from 
Table 2. This property is consistent with our numerical re- 
sults as shown in Figures 1 and 2. 

In Figures 3 and 4, typical distributions of velocity per- 
turbations Sv r and Svg are displayed for the r-modes of 
slowly and rapidly rotating stars, respectively. Here we de- 
fine, 1 



8v r = 5u r 



Svg = r&u 



(17) 
(18) 



Here, after discarding data near the mass-shedding limit, we 



The equilibrium models shown in the figures are taken 
from the sequence b in Table 1. The rotation parameters of 
the stars shown in Figures 3 and 4 are chosen as r p /r e = 0.96 
(T/\W\ = 0.01248) and r p /r e = 0.70 (T/\W\ = 0.1035), re- 
spectively. In each figure, the amplitudes of the eigenfunc- 
tions are shown against the surface-fitted radial coordinate 
r, (0 < r* < 1) for five different constant values of #*, whose 
values are given in the figures. Figures 3 and 4 show that 
the relativistic r-mode oscillations of rotating stars are very 
similar to the corresponding Newtonian modes. Moreover, 
comparing Figure 3 with Figure 4, we observe that the ba- 
sic properties of the r-mode oscillations of relativistic stars 
depend only weakly on the magnitude of the rotation. How- 
ever, it is important to emphasize the following properties; 
the horizontal motion of the fluid due to the r-mode os- 
cillations is more concentrated near the surface of the star 
for a rapidly rotating star than for a slowly rotating star. 
Similar behaviors are observed in the r-mode oscillations of 
Newtonian stars (Karino et al. 2000). 

For slowly rotating stars, the r-mode oscillations of 
Newtonian stars do not excite the fluid motion in the ra- 
dial direction, while in the relativistic star, the fluid element 

1 Note that these are not exactly the components of 3 velocity on 
the orthonormal frame which is naturally defined on our space- 
time. 
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can move in this direction due to the effect of general rela- 
tivity (Lockitch, Andersson & Friedman 2000). This prop- 
erty can be confirmed in Figure 3, which shows that the 
perturbations of the radial velocity Sv r for a slowly rotat- 
ing relativistic star have non-vanishing amplitudes. On the 
other hand, as already shown in Newtonian studies (Provost, 
Bcrthomieu & Rocca 1981; Saio 1982), the amplitude of 5v r 
for the r-mode oscillations becomes large as a rotation pa- 
rameter T/|W| of the equilibrium star is increased. This 
fact is also observed in Figure 4. In summary, our results for 
modal properties are consistent with previous results that 
properties of pure-axial perturbations for the r-mode oscil- 
lations cease to appear for a highly relativistic star and/or 
for a rapidly rotating star. 



4 DISCUSSION AND SUMMARY 

4.1 Discussion 

It is the most crucial thing to know how far we could rely 
on the results obtained by using the Cowling approximation. 
Since no one has succeeded in developing a formulation for 
the full analysis of r-mode oscillations of rapidly rotating 
stars in general relativity, what we can do is to estimate the 
accuracy by considering the results obtained for the New- 
tonian configurations and those for slowly rotating general 
relativistic stars. 

The frequencies of m = 2 r-mode oscillations for the 
Newtonian polytropes are shown in Table 3. In this table, 
the ratios of the r-mode oscillation frequency to the angular 
velocity in the full problem and in the Cowling approxima- 
tion are tabulated both for slowly rotating polytropes and 
for rapidly rotating polytropes with the polytropic indices 
N = 0.5 and 1. As seen from this table, the values of the 
Cowling approximation are to within a several percent from 
the exact values. Since there occurs no density perturbation 
in the r-mode oscillation of spherical Newtonian stars, the 
Cowling approximation gives very accurate values for slowly 
rotating Newtonian configurations irrespective of the poly- 
tropic indices. This kind of behavior obtained by using the 
Newtonian Cowling approximation is not directly applied 
to the general relativistic configurations because the basic 
equations in general relativistic oscillations are not the same 
as those in the Newtonian models. However, it should be 
noted that the results of the scaled eigenfrequencies a /Q. 
obtained by the Cowling approximation are located within 
a several percent from those of the full problem. Together 
with the errors for general relativistic spherical stars dis- 
cussed in the previous section, we may as well believe that 
the results obtained in this paper for the r-mode oscillations 
of rapidly rotating general relativistic polytropes would be 
also well approximated ones. 

In general, the master equations for oscillations of a 
stationary and axisymmetric rotating star can be reduced 
to a single two-dimensional second-order partial differential 
equation as shown below. However, the type of the reduced 
second order partial differential equation cannot be known 
beforehand because it depends on the oscillation frequency 
of the mode (see, e.g., Balbinski 1985; Skinner & Lindblom 
1996; Lindblom 1997). In other words, it implies that the 
type of the equation cannot be determined before some so- 
lutions of the modes are obtained. This is an important fact 



Table 3. .R-mode frequencies of uniformly rotating Newtonian 
polytropes. The results for the full problem and those obtained 
by the Cowling approximation are tabulated for N = 0.5 and 1 
polytropes. Here r p /r e is the axis ratio of the rotating configu- 
rations, where r p and r e are the polar radius and the equatorial 
radius, respectively 



Method 


N 


r p /r e 


a/Q 


Full 


0.5 


0.98 


1.320 


Cowling 


0.5 


0.98 


1.319 


Full 


0.5 


0.67 


1.099 


Cowling 


0.5 


0.67 


1.050 


Full 


1.0 


0.98 


1.323 


Cowling 


1.0 


0.98 


1.321 


Full 


1.0 


0.67 


1.159 


Cowling 


1.0 


0.67 


1.116 



because the properness of the boundary value problem for the 
partial differential equations is deeply related to the type 
of the differential equation. If the equation is of the ellip- 
tic type, the boundary value problem is a properly posed 
one. However, if it is hyperbolic, one might need to devise 
a proper treatment of the boundary value problem for the 
hyperbolic equations. 

In our relativistic Cowling problem, we can derive the 
following second order partial differential equation in terms 
of the pressure perturbation and consider its type. From the 
basic equations (7) and (11)-(13), we can obtain 



,d 2 8p ^d 2 8p ^d 2 8p „,d8p dSp . 



= , (19) 



A 



(a - raQ) 2 
e+p u l D(o — mQ,) 2 1 £7 — u> + ujv 



B = 



1 dlaFdhiH 



(20) 



e+p v? D(a ■ 
f dlnFdlnH 
I dr 86 4 



- mQ) 2 r 2 Q. — lj 
d\nF d\nH \ 
~d~6 dr~ J ' 



(21) 



C = 



D 



e + p 


'l/ f 


D(a 


-mQ.) 2 


r 2 




f (- 


— m 


nf 


v 2 


dlnFdlnH\ 


\Q- 


UJ + 


- OJV 2 


Q — UJ 


dr 


dr j ' 


1 


I 


' (<T- 

n- 


-mQ) 2 

UJ + OJV 2 






a — mfl \ 




2 2v 

v e 


-la 


(dlnF dlnH 


1 


dlnFdlnH 



(22) 



O, — uj 



dr dr 



(23) 



where / is a function of the first order derivative of Sp, Sp 
itself, and coordinates (r, 9) whose form is not essential for 
the following discussion (see, also, Ipser & Lindblom 1992). 
The type of the above partial differential equation (21) 
at each mesh point can be determined by the following quan- 
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tity Q: 

Q = B 2 - 4AC 



(24) 



For a uniformly rotating Newtonian barotropic star, the 
quantity Q is given by 

16 



Q = 



l 



r 2 AQ 2 - (a- - mQ) 2 



(25) 



Accordingly, the master equation in the Newtonian limit is 
hyperbolic (elliptic) if the frequency satisfies 4Q 2 > (a — 
mQ,) 2 (AQ 2 < (cr — mQ) 2 ). It is important to note that oscil- 
lation frequencies of all inertial modes of uniformly rotating 
Newtonian stars that have been obtained so far are in the 
range of —20. < a — mQ < 2fi, implying that Eg. 11911 is the 
hyperbolic (cf. Greenspan 1969). 

As for all general relativistic r-mode oscillations ob- 
tained in this paper, values of the quantity Q are also posi- 
tive at every point inside the star. It follows that the equa- 
tion is again hyperbolic. Although there is no mathematical 
proof, therefore, it seems that solutions of the inertial modes 
are basically given as solutions of partial differential equa- 
tions of the hyperbolic type. Note that the situation would 
be more complicated if non-barotropic stars were consid- 
ered. In non-barotropic stars, the inertial modes behave like 
g-modes which are solutions of an elliptic equation when the 
rotational velocities of the stars are sufficiently slow (see, 
e.g., Yoshida & Lee 200b). 

4.2 Summary 

In this paper we have developed a numerical scheme for ob- 
taining the r-mode oscillations of rapidly rotating relativis- 
tic stars. In the present scheme, we neglect all the metric per- 
turbations and only take account of the dynamics of the fluid 
motion in the fixed background spacetime of the star (the 
relativistic Cowling approximation) . We also assume stars to 
be barotropic under the assumption of adiabatic oscillations. 
Our numerical scheme is based on the Yoshida-Eriguchi for- 
mulation (Yoshida & Eriguchi 1997) and is a general rel- 
ativistic extension of the Karino- Yoshida- Yoshida-Eriguchi 
numerical scheme for determining oscillations of rapidly ro- 
tating Newtonian stars (Karino et al. 2000). With this new 
numerical scheme, the frequencies of the r-mode oscillations 
are obtained as functions of the ratio of the rotational en- 
ergy to the absolute value of the gravitational energy T/\W\ 
along sequences of polytropic equilibrium stars whose ratios 
of the pressure to the total energy density at the center of 
the stars and the polytropic index are kept constant. It is 
found that scaled oscillation frequencies a/Q are decreasing 
functions of T/II4 7 ! and are given as nearly linear functions 
of T/|W|. As already found in studies on the relativistic r- 
mode oscillations of slowly rotating stars, we observe that 
the r-mode oscillation frequencies are sensitive to change of 
the relativistic factor of the stars but not to the change of 
stiffness of the equation of state. The present results are con- 
sistent with the previous results of the r-mode oscillations 
both for rapidly rotating Newtonian stars and for slowly 
rotating relativistic stars. 

It should be noted that three simplifications about the 
basic properties of the relativistic r-mode oscillations of 
rapidly rotating general relativistic stars have been em- 
ployed in the present investigation; namely, the use of (1) 



the relativistic Cowling approximation, (2) the rigid rota- 
tion law, and (3) the barotropic equation of state. Although 
it is known that the Cowling approximation can produce ap- 
proximate eigenfrequencies with acceptable errors and basic 
modal properties similar to the exact ones for such oscil- 
lation modes as the fluid motion is dominating, the exact 
treatment of full equations for obtaining oscillation modes 
of rapidly rotating relativistic stars is required for under- 
standing accurate solutions for the oscillation modes and 
for exploring purely general relativistic properties such as 
damping times of oscillations due to gravitational radiation 
reactions. A great challenging problem to find a method for 
obtaining and solving exact pulsation equations for rapidly 
rotating relativistic stars has been left unsolved yet. 

An extension of the present method to differentially ro- 
tating stars and/or to non-barotropic stars would be worth- 
while in the following context: nascent neutron stars are ex- 
pected to rotate differentially and the effect of differential 
rotation could strongly affect the frequencies and the modal 
properties of the r-mode oscillations (for Newtonian models, 
see e.g. Karino et al. 2001). In order to have a correct sce- 
nario of spin evolutions of young neutron stars due to the r- 
mode instability, therefore, we have to know how the r-mode 
oscillations would be changed under the effect of differential 
rotation. As for the r-mode oscillations of non-barotropic 
stars, as mentioned in Section I, we do not yet know whether 
the r-mode oscillations similar to those of Newtonian stars, 
whose eigenfrequencies are discrete and eigenfunctions are 
regular, can exist in relativistic non-barotropic stars. Fu- 
ture investigations on the relativistic r-mode oscillations of 
rapidly rotating non-barotropic stars in the Cowling approx- 
imation might give us useful clues for understanding of the 
full r-mode oscillations of non-barotropic relativistic stars. 
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Figure 1. Scaled frequencies a/Q. of the m = 2 r-mode oscil- 
lation of rotating stars with the polytropic index TV = 0.5 (se- 
quence a and sequence b) are plotted as functions of the ratio of 
the rotational energy to the absolute value of the gravitational 
energy, T/|W|. The frequencies of the r-mode oscillations for the 
sequence a and the sequence b are shown by dashed and solid 
curve, respectively. 
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Figure 2. Same as Figure 1 but for sequences c (dashed curve) 
and d (solid curve). 
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Figure 3. Eigenfunctions 8v r (top) and Svq (bottom) of the 
m = 2 r-mode oscillation for a slowly rotating star. The equi- 
librium configuration belongs to the sequence b and its rotation 
parameter is taken to be r p /r e = 0.96 (T/\W\ = 0.01248). Distri- 
butions of the eigenfunctions for several values of 8* are displayed 
as functions of the surface fitted radial coordinate r„ . Normaliza- 
tion of the eigenfunctions is taken so as to be Svg = 1 at the 
stellar surface on the equator. Note that 5v r for 9 t = 7r/10 and 
it/2 are very small and cannot be seen in this figure. 



Figure 4. Same as Figure 3 but for a rapidly rotating star whose 
rotation parameter is taken to be r p /r e = 0.70 (T/114 7 = 0.1035). 



